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The concordance model of cosmology contains several unknown components such as dark matter 
and dark energy. Many proposals have been made to describe them by choosing an appropriate 
potential for a scalar field. We study four models in the realm of loop quantum cosmology: the 
Chaplygin gas, an inflationary and radiationlike potential, quintessence and an anti-Chaplygin gas. 
For the latter we show that all trajectories start and end with a type II singularity and, depending 
on the initial value, may go through a bounce. On the other hand the evolution under the influence 
of the first three scalar fields behaves classically at times far away from the big bang singularity and 
bounces as the energy density approaches the critical density. 
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PACS numbers: 04.60.Pp,98.80.Cq 

I. INTRODUCTION 

It is generally believed that our Universe started with 
an inflationary phase followed by a radiation and mat- 
ter dominated era. However, classical cosmology is not 
able to tackle the problem of the initial conditions of 
the universe. One of the possible solutions to this prob- 
lem is that our expanding Universe was preceeded by 
a contracting phase. But powerful singularity theorems 
based on classical general relativity (GR) forbid such a 
behavior unless one assumes a form of matter that vio- 
lates the positive energy conditions or modified versions 
of gravity (see e.g. [l|, Q)- On the other hand it is be- 
lieved that quantum gravity should solve this problem 
by generating ideal conditions for the genesis of our uni- 
verse. Several proposals such as the pre-big bang string 
cosmologyQ and the ekpyrotic/cyclic models[l,[3l mod- 
ify dynamics with (perturbative) quantum gravitational 
effects but have so far had limited viability. 

A generic nonsingular evolution through the big bang 
can only be achieved if nonperturbative effects of quan- 
tum gravity are incorporated. One of the leading non- 
perturbative background independent approach is loop 
quantum gravity (LQG)0-[1]. One of the main predic- 
tions of LQG is that the underlying geometry is discrete. 
The application of the quantization methods of LQG to 
homogeneous spacetimes results in what is known as loop 
quantum cosmology (LQG) [§-15|. The results of LQG 
not only provide new insights into the quantum struc- 
ture of spacetime near the big bang singularity but also 
remove this singularity by extending the time evolution 
to negative times. It has been rigorously shown 16- isj 
that the evolution of contracting semiclassical universes 



passes through the quantum regime and expands to semi- 
classical universes. This nonsingular bounce stems from 
the fact that the dynamics in LQG is governed by a dis- 
crete quantum difference equation in qiiantu m g eometry. 
On the other, it can be shown [3, E3, 0, Eo] that an 
effective Hamiltonian on a continuum spacetime can be 
found which approximates well the quantum dynamics 
(for a Wheeler-DeWitt analog see 2l|). The modifica- 
tion arising from nonperturbative effects to the classical 
Friedmann equation includes a term p^/pcrit, where p is 
the energy density and pcrit denotes the critical density 
of the order of magnitude of the Planck density. Since 
this term is negative the evolution bounces whenever the 
energy density reaches a density close to the Planck den- 
sity. 

The viability of the bounces for more general mat- 
where it 
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ter sources has been studied in e.g. 
was shown that the behavior of solutions with inflation- 
ary and negative potentials are nonsingular, respectively, 
where solutions of exponential potentials are analyzed. 
Moreover it was shown that for negative potentials the 
inner boundary also appears, corresponding to the clas- 
sical rccoUapse, which leads to solutions having cyclic 
behavior. In 23, 13| the authors studied the role of LQG 
effects in the Ekpyrotic/Gyclic model in Bianchi type I 
models and showed that the universe undergoes multiple 
small bounces and the anisotropic shear remains bounded 
throughout the evolution. 

In this work we are interested in potentials which play 
a major role in classical cosmology. We first introduce 
the effective dynamics in LQG in Sec. |lll In Sec. IIIII we 
give a short overview of conditions for singularities oc- 
curing in Friedmann-Robertson- Walker (FRW) cosmolo- 
gies. In Sec.|lV]we study the Ghaplygin gas and in Sec.lVl 
we study the robustness of the bounce for a scalar field 
which has the property of being inflationary at small 
times and radiationlike at later times 28143(1. Sec. IVII 



is devoted to the anti-Ghaplygin gas and Sec. I VIII to a 
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quintessence model which models dark energy. We con- 
clude with Sec. IVml 



II. EFFECTIVE DYNAMICS IN LQC 

LQG is a canonical quantization of gravity based upon 
Ashtekar connection variables. The phase space of classi- 
cal GR in LQG is spanned by a SU(2) connection and 
a densitized triad Ef on a 3- manifold M , which are two 
conjugate variables encoding curvature and spatial ge- 
ometry, respectively. Likewise, LQC is a canonical quan- 
tization of homogenous spacetimes such that the phase 
space structure is simplified, i.e., the connection is de- 
termined by a single quantity labeled c and likewise the 
triad is determined by a parameter p. For the spatially 
flat model of cosmology, the new variables are related to 
the metric components of the (FRW) universe through 
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where 7 is the Barbero-Immirzi parameter which is set to 
be 7 sa 0.2375 by the black hole entropy considerations 
(3H . Classically in terms of the connection-triad variables 
the Hamiltonian constraint is given by 
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with K = SttG (where G is Newton's gravitational con- 
stant) and "Hm being the matter Hamiltonian. The com- 
plete equations of motion are derived from Hamilton's 
equations x = {x,Hci} for any phase space variable x, 
and by enforcing that Hci should vanish. The variables 
c and p are canonically conjugate with Poisson bracket 
{c,p} = 7^/3. The classical Friedmann equations are 
obtained through a substitution of these relations into 
the Hamiltonian constraint ©. 

The basic variables of LQC are the component of the 
densitized triad and the holonomies along straight edges: 
hi{n) = exp(//cTi) = cos(/ic/2) + 2 sin(/ic/2) t^, where 
is related to the Pauli spin matrices through r,; = —iai/2. 
The dimensionsless parameter fi represents the physical 
length of the edge and is arbitrary. As such, the almost 
periodic functions exp(i/ic/2) =: M^{c), G R, can be 
chosen to be the elementary variables of LQC. The op- 
erator p corresponding to the component of the densi- 
tized triad acts by differentiation and is diagonalized by 
Mii{c). In a canonical setting, the dynamics is imple- 
mented completely by the Hamiltonian constraint. Upon 
quantization, the Hamiltonian constraint is promoted to 
an operator using Thiemann's trick [3^ . 

It has been shown that the underlying dynamics in 
LQC is governed by a discrete difference equation in 
eigenvalues Vp, of the volume operator V in quantum 



geometry (see e.g. Q). However, an effective Hamil- 
tonian description on an continuum spacetime can be 
constructed using semiclassical states, which has been 
shown to very well approximate the quantum dynamics 
[igI \v\ . The analysis of the quantum Hamiltonian us- 
ing semiclassical states belonging to the physical Hilbert 
space reveals that a backward evolution of our expand- 
ing phase of the universe leads to a bouncing solution 
into a contracting branch 33[ . The expectation values of 
the Dirac observables allows us to investigate to quantify 
the difference between the quantum and classical dynam- 
ics. It turns out that quantum geometric effects become 
dominant only when the energy density p of the universe 
is of the order of the critical density Pcnt [17|, l3J] and 
classical general relativity is a good approximation to 
the quantum dynamics when p <^ Pmt- The effective 
equations for the modified Friedmann dynamics can be 
derived from the effective Hamiltonian constraint with 
loop quantum modifications. The effective Hamiltonian 
constraint, to leading order, is given by [19 1 
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Heff = —asiT?{pc) -{- Hm ■ (3) 



where fi = ^J3V^/2\p\ [l7|. 

In this work we will be mainly interested in the matter 
Hamiltonians corresponding to a massive scalar field (f> 
with momentum H^ and potential V{(j)): 



nu = l^+p'/'vi^). 



(4) 



The energy density and pressure of the scalar field are 
given by 

P = P4> = l^^ + V{cj,), p^ = l^^-V{cj,) . (5) 

The scalar field satisfies the stress-energy conservation 
law: 



p^ + 3-(p0 = 



(6) 



The modified Friedmann equation for p is obtained 
with Hamilton's equations 

, , JK d 2a _ X 

p = {p, rtcff I = — Tr^rtcff ^ — am{p.c) cos(/ic) (7) 
3 oc JIJ. 

which on using Eq. ([T]) implies that the rate of change of 
the scale factor is given by 

I 



a = sm{fic) cos(/ic) 
7// 



(8) 



Furthermore, the vanishing of the Hamiltonian constraint 
implies 
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Combining Eqs.® and (|9]) yields the effective Fricdmann 
equation for the Hubble rate H = a/a 
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3 V Pcrit 

with the critical density given by 
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where ppi = l/(fi,G^) is the Planck density. The modified 
Fricdmann equations provides an effective description for 
LQC which approximates the underlying quantum dy- 
namics very well. The ^^-modification of the Fricdmann 
equation arises from nonperturbative quantum geometric 
effects. Since the modified term is negative definite, the 
Hubble parameter vanishes when p = pcrit and the uni- 
verse experiences a turn-around in the scale factor. For 
p ^ Pcrit, the modifications become negligible such that 
the standard Fricdmann equations are recovered. In ad- 
dition, it should be noted that pcrit is of quantum origin 
since Ppi oc 

The modification arising in Eq. ([TU]) constrains the 
Hubble parameter and the energy density to be bounded 
from below and above: 



KPcrit 
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] and p e [0, Pcrit] • (12) 



Moreover, using the conservation law ([S]), its time 
derivative can be cast into 



if = -^(P + P0) (l- — 

^ \ Pcrit 



(13) 



The Hubble parameter can be expressed in terms of the 
scalar field such that 
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where expansion corresponds to the upper sign and con- 
traction to the lower sign. 

Apart from the case of a free scalar field it is difficult 
to find analytical solutions. For this reason we will draw 
phase portraits showing the qualitative behavior of the 
numerical solutions. Using the Hamiltonian ^ the four- 
dimensional parameter space (c(t),p(t), 0(i), </'(i)) can be 
reduced by one unit by expressing one of these variables 
by the other three. Following [2^ we will display a phase 
portrait consisting of the variables cj) and </>. The quantum 
turn-arounds will be represented as solid lines showing a 
boundary for the solutions. Once a trajectory reaches 
such a boundary the sign of the Hubble rate changes, 
thus indicating a turnaround. Since our potentials are 
complicated we will also show the phase portraits illus- 
trating the time evolution of the kinetic and potential 
energy. 



III. COSMOLOGICAL SINGULARITIES 

Over the last few years, the zoo of cosmological singu- 
larities has become considerably more extensive. Beside 
the traditional singularities known as big bang an big 
crunch, there also exist the big rips and sudden singu- 
larities. In this section we provide a catalog of relevant 
singularities in FRW cosmologies 35-3l|. All singulari- 



ties are classified by means of both the kinematic (scale 
factor a) and the dynamical (energy density p and pres- 
sure p) behavior . 

Big bang and big crunch: the most basic of the cos- 
mological milestones are big bangs and big crunches, for 
which the scale factor a{t) — ^ at some finite time as 
we move to the past or future. Also both the energy 
density p and curvature invariants diverge. Dynamically, 
the Null Energy Condition (NEC hereafter), p + p>Q, 
is always satisfied. 

Big rip or type I singularity: a big rip is said to occur 



40[. This is accom- 



if a{t) — > 00 at some finite time [SE 
panied with a divergence of the energy density, pressure 
and curvature invariants. These singularities always vi- 
olate NEC and all other energy conditions such as the 
WEC (NEC & p > 0), SEC (NEC & p + 3p > 0) and 
DEC (p±p > 0) [37|. The type I singularity emerges for 
phantom-like equations of state: w < —1. 
Sudden or type II singularity: this extreme event is char- 
acterized by a finite value of the energy density but an 
associated divergence of pressure at finite value of the 
scale factor. Due to the divergence of pressure, the Ricci 
curvature scaler R diverges. 

Type III singularity: like type II singularities, but the 
energy density and pressure diverge, causing a blow up 
of curvature invariants at some finite time. The type HI 
appears for the model with p 
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Type IV singularity: higher-order time derivative of the 
scale factor a diverge at finite time, while the scale fac- 
tor itself remains finite. None of the energy density or 
pressure blows up in this case. The type IV singularity 
appears in models characterized b y p = —p — f{p), where 
/(p) can be an arbitrary function 35 [. 



IV. CHAPLYGIN GAS 



The Chaplygin gas was introduced to cosmology in [43 1 
to describe the transition from a universe filled with dust- 
like matter to an exponentially expanding universe. This 
gas IS a perfect fluid which has the following equation of 
state: 



A 
p 



(15) 
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where ^ is a positive constant. Energy conservation re- 
quires that 



where B is an integration constant. For positive B and 
small a we get a universe dominated by dust-like matter: 



B 

T 



For large a it turns out that the universe is of the de Sitter 
type with a cosmological constant \/A: 



P' 



a" > 
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'a' 



The potential corresponding to this equation of state is 
given by [43 



Vch(0) 
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= 2"^/^ f cosh\/3K0 
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A generalization of this gas has been introduced in [4J], 
where the equation of state is given by 

A 

P = — ^' 
P 



where a is a positive constant. The requirement that 
the sound velocity not exceed the speed of light yields to 
the bound < a < 1. The potential for the scalar field 
corresponding to this equation of state reads 



^gch(0) - Vo cosh(^/;^/30) + cosh(V^/3(/.) ' 



where /3 = y^3{a + l)/2. Leaving both A and a free, 
the latest cosmological and astrophysical constrain these 
parameters to the following la confidence level [iHl 

a = -Omto.ll and As = 0.73l^:o^, 

where Ag = Aj i^A-VB) and B is an integration constant. 
However, we checked numerically the influence of a and 
it turns out that different values of a do not change the 
behavior of the bounce. The reason is that, since a ^ 
implies |0| — 00, we have 

ltcH'-K)Cxp(^/;^/3|0|)^. 



Fig. [T] shows the phase portrait of the variables 
{(j){t),(j>{t)) for four different initial values. All trajec- 
tories start from the point (0,0) for t — > —00 where the 
energy density vanishes. Also, the Hubble rate is nega- 
tive but close to zero (see Fig. [5]). As can be seen from 
Fig. [3] the potential energy is the dominant contribution 
for times far away from the bounce. As the energy den- 
sity starts to grow the trajectories depart from (0, 0). At 
time t = —4.3 the kinetic enery vanishes and the poten- 
tial energy reaches a local maximum, which can also be 
seen from the plateau in the Hubble rate. The evolu- 
tion then reaches the point at the boundary where the 
universe bounces. The energy density is highest at this 
point (denoted by the dot in Fig. [3]) and the dominant 
contribution comes from the kinetic energy. Moreover, 
as a manifestation of the bounce the sign of the Hubble 
rate changes and becomes positive. Then the evolution 
reaches a second plateau at time t = 1.9 where the kinetic 
energy vanishes and the potential has its global maxi- 
mum. For i — > 00 all trajectories go to the point (0, 0) 
of the phase portrait and the Hubble rate decreases with 
the same rate as radiation, i.e. H ^ {2t)~^. 




FIG. 1: Phase portrait for the Chaplygin gas with poten- 
tial (|16p for y/A = 10~^. The thick line shows the boundary 
indicating turn-arounds, the thin lines show the solutions of 
Eq. (|14p for different initial data. 



V. INFLATIONARY AND RADIATIONLIKE 
POTENTIAL 



Thus, different values in the parameter space (A, a) only 
change the magnitude of the potential and not its shape. 
From now on we only consider the case a = 1. 



In this section we study a scalar field with a scalar 
filed potential which can be seen as a modification to 
the Chaplygin gas. The general motivation behind this 
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FIG. 2: Graphs of the Hubble rate and its time derivative 
for the Chaplygin gas as a function of time. The dashed hne 
represents tlie classical solution and the solid lines the solution 
from LQC. The bounce occurs at t = 0.52, where the Hubble 
rate changes sign. 




FIG. 3: Phase portrait of the kinetic energy -Ekin ~ 4>'^ /"^ 
and the potential energy Vch(<;/>). The dot corresponds to the 
bounce. The kinetic energy vanishes at t = —4.3 and t — 1.9. 



modification is the fact that the potential can be used 
to model both radiation and inflation (see also 4^ 47|). 
Let us consider the energy density [28|-|30|: 



A 



B 



a4(i+«) 



1 + a < 0, 



(17) 



where A, B and a are constants. For early times the 
energy density is inflationary: 

and for late times radiationlike: 

p^\, B/a^^^+"^ :^ A. 

Such a behavior can be modeled by a scalar field with 
the following potential 



Vir(0) = 



Vo 



coshi+= (-fc(l + «)(/)) 

+ 2 cosh^TT^ (-fc(l + a)(p) 



(18) 



This potential shares many similarities with the Chaply- 
gin potential (fTH]). On the other hand, while the potential 
energy of the Chaplygin gas is the dominant contribution 
to the energy density at late times and just after (respec- 
tively before) the bounce (see Fig. [3]), Vm{4>) is always 
at least one order of magnitude smaller than the kinetic 
term. So, instead of going to zero as can be seen in Fig.[l] 
(f) ±oo as i — )■ ±oo (see Fig. S]). Moreover, the Hub- 
ble parameter H decreases at a radiationlike rate, i.e. 
H ^ 1/ {2t) . As the evolution approaches the classical 
singularity both the kinetic and potential energy densi- 
ties increase and approach the critical density. This is 
when LQC modifications come into play such that the 
evolution goes through a bounce. This point is reached 
when the evolutions represented by the thin lines in Fig. 2] 
touch the boundary shown as thick lines. 




FIG. 4: Phase portrait for a scalar field with potential (|18p for 
Vo = 10"^. The thick lines represents the boundaries where 
the solution of Eq. (|14|) undergo a bounce. 



VI. ANTI-CHAPLYGIN GAS 

The anti-Chaplroin gas was introduced in the context 
of cosmology in [48|. The pecularity of this gas is that 
its equation of state is given by 

A 

P = — , 
P 

which can be modeled by a scalar field with potential 
VachW = Vo UnhiV3^\cj)\) - smh-\V3^\(b\)) . (19) 



In a cosmological context such a gas leads to a big brake 
singularity caused by the divergence of of higher deriva- 
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tives of the Hubble rate. This singularity occurs at a 
finite value of the scale factor where the Hubble rate van- 
ishes. Since the second time derivative of the scale factor 
diverges the Ricci scalar R also diverges. Moreover, while 
the energy density remains finite the pressure diverges. 
Such singularities are called sudden of type H singularity 

The energy density is near zero when the type H sin- 
gularity occurs such that the modifications arising from 
LQC are not able to avoid this divergence. This can 
be seen in Fig. [5] where every solution converges toward 
(j) = 0. The reason why the energy density does not 
diverge is because the kinetic term in Eq. ^ cancels 
the divergence from the potential. However, the pres- 
sure diverges because the potential is unbounded from 
below. On the other hand LQC is able to remove the Big 
Bang singularity occuring at times represented by dots in 
Fig. [5] From a backward evolution perspective, not ev- 
ery solution goes through a bounce because, depending 
on the initial value, the third singularity is reached. This 
singularity is also of type H because the energy density 
is finite but the pressure diverges. As before, LQC is not 
able to resolve it and the evolution stops. 

In sum, the evolution of a universe filled with an anti- 
Chaplygin gas starts and stops at a type H singularity 
when the point (/) = is reached. Depending on the 
initial value, it may go through a bounce. 




FIG. 5: Phase portrait for the anti-Chaplygin gas with po- 
tential (|19p for Vb = 10"*. The thick line represents the 
singularity Vach(0) — >■ — oo, the thine lines the solution of 
Eq. (|14|) and the dots the bounce. The evolution starts and 
ends at a type II singularity. 



VII. QUINTESSENCE 

Recent observations of the anisotropy of the cosmic 
mircowave background (CMB) [i^ together with the 
power spectrum of the large scale structure (LSS) [H^ 
and the magnitude-redshift relation of the supernovae la 
51 , l^l all indicate that the current mean energy density 
Ptot of the Universe consists not only of radiation, bary- 
onic and dark matter, but also of a dominant component 
of negative pressure form which is called dark energy. An 
explanation for the missing energy is quintessence where 
the dark energy density is identified with the energy den- 
sity (associated with a negative pressure p^) arising 
from a scalar (quintessence) field 4>. It is possible to con- 
struct scalar field potentials V which lead to a constant 
equation of state = p^/p^. The form of such poten- 
tials depends on the energy components of the model. 
Potentials with three components (radiation, matter and 
qiuntessence) can only be modeled for special values of 
uu. The exact quintessence potential for : 



I reads 



Vq(0) - 



Vo 



[r]smh{B(l)) + cosh{B(p) - 1]' 



(20) 



where the potential strength Vq, respectively B and rj are 
given by 



Vo = 



8 r^^o^ 

3 "02" 



-Po, 



B = 




77 = 2 



^(21) 

and = 1 — Oiad — Om, where we utilize dimensions- 
less density parameters fti = Pi/ Po with po = /SttG. 
We use a model consistent with Wilkinson Microwave 
Anisotropy Probe (WMAP) 5-year data [i^. The poten- 
tial ([20| and therefore the cosmic evolution is governed 
by two very different energy scales: the huge Planck mass 
rripi and the much smaller energy density po- Explicitly, 
one derives from (pil)) for (f) ^ 



1 



and, respectively, for 



oo 



FQ((/))^exp(-2i?^) , 



(22) 



(23) 



which is in accordance with the general behavior of a 
quintessence potential. 

Numerical solutions for the time evolution in LQC with 
a quintessence potential are shown in Fig. [71 respectively. 
Fig. [51 and the phase potrait, consisting {(f), (j)), for differ- 
ent initial values is presented in Fig.[6l As in the previous 
cases both the Hubble parameter H and the energy den- 
sity p are bounded and subject to the constraints (|12l) . 
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As can be seen from Fig. [7] the Hubble rate starts with a 
small negative value. For times far away from the bounce 
the potential energy is the dominant contribution, as can 
be seen in Fig. |S1 All trajectories start from 4> — > ±00 for 
t — >■ — cxD, cf. Fig. ini When the kinetic energy reaches a 
value such that the energy density becomes comparable 
to pcrit the magnitude of the Hubble rate starts increas- 
ing and quickly becomes zero at p = pcrit or equivalently 
aX t = —0.86. As shown in Fig. [5] the evolution than 
reaches a point at the outer boundary, where the universe 
bounces. The bounce implies a change of the sign of the 
Hubble rate, cf. Fig. [71 Immediately after the bounce the 
universe expands quickly and reaches a plateau at time 
t = —0.5, where the potential energy reaches a global 
maximum, whereas the kinetic energy (j)/2 vanishes. As 
presented in Fig. [3 for t — >■ 00 the Hubble parameter 
decreases at a radiationlike rate, i.e. H ^ {2t)^^ and all 
trajectories starting from — > ~oo for t — > ~oo go back 
to (/) — ^ —00 and, respectively, all trajecories coming from 
— >■ 00 end in — > 00 for large times. Thus, there exist 
two independent sectors in the phase diagram, cf. Fig. [51 













f 1 















FIG. 6: Phase portrait with the quintessence potential (|20|l . 
The outer boundary (thick line) corresponds to p = /Ocrit, the 
thin lines show the solutions of p4|l for different initial data. 



VIII. DISCUSSION AND CONLUSIONS 

As an attempt to solve the shortcomings of the concor- 
dance model of cosmology several models of scalar fields 
have been proposed which interpolate between two stages 
of the evolution of our universe. We studied the influence 
of three types of scalar fields of cosmological interest. 




-10 -5 5 10 



0.2 
0.0 
-0.2 
-0.4 
-0.6 
-0.8 
-1.0 




-10 -5 5 10 



FIG. 7; Behavior of H and H with the quintessence model 
potential H20p . The dashed line shows the classical solution, 
whereas the solid lines show the solutions obtained from LQC. 
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FIG. 8: Phase portrait of the kinetic energy 0^/2 and the 
potential energy VQ((/i), where the dot corresponds to the 
bounce. The kinetic energy vanishes at t = 0.5. 



nameley the Chaplygin gas, a modificated version of it 
and quintessence, and one more exotic type called the 
anti-Chaplygin gas. While the first type models a uni- 
fication of dark matter and dark energy, the second one 
interpolates between an early inflationary phase and ra- 
diation. Contrary to quintessence which was introduced 
as an effort to describe dark energy in terms of a scalar 
field, the anti-Chaplygin can be considered as a toy model 
without any direct application to cosmology. 

We presented the solutions to the Friedmann equations 
in LQC for these four models. We showed that the evo- 
lution of the first three models (Chaplygin gas, modi- 
fied Chaplygin gas and quintessence) follows the classical 
path until it approaches the critical density, where the 
modification to the Friemann equation gains in impor- 
tance. As this modification is negative definite the evo- 
lution bounces and the Hubble rate changes sign. Some 
time after the bounce the evolution follows once again 
the classical trajectory. We showed that, while all the 
origin in the (0, (/))-phase diagram acts as an attractor 
for the Chaplygin gas, the solutions of the modified ver- 
sion converge toward — and cf) — > ±00. Quintessence 
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behaves in a similar way, except that there are two in- 
dependent sectors in the (0, (/))-phase diagram such that 
trajectories with positive respectively negative initial (p 
always stay positive respectively negative. The situation 
is radiacally different for the anti-Chaplygin gas where 
every trajectory starts and ends with a Type II singu- 
larity. Depending on the initial data the evolution may 
go through a bounce, however LQC is, as expected, not 
able to remove these Type II singularities. Because of 



this very fact there are also two independent sectors for 
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